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^^ ; Abstract 

r~| ■ We first show that the rotatmg mass matrix hypothesis suggested ear- 

Q4 her, where the massive eigenvector of a rank-one mass matrix changes 

o ' with renormalization scaie, is consistent with the latest experimental 

^ • data on fermion mass hierarchy and mixing, including the CP violat- 

r^i ' ing KM phase. We obtain thereby a smooth trajectory for the massive 

eigenvector as a function of the scale. Using this trajectory we next 

^ ■ study Higgs decay and find suppression oi T{H — )• cc) compared to 

VO . the standard model predictions for a range of Higgs masses. We also 

J^ I give limits for flavour violating decays, including a relatively large 

^vq ■ branching ratio for the r^^u"*" mode. 

O! 

O- 1 Introduction 

_ ^ I Even though the standard model is a tremendously successful theory, there 

^ ■ are many aspects of it which are not understood. Among the most puzzling of 

these is the hierarchy of masses across generations. The rotating mass matrix 
picture purports to give a natural answer to the mass hierarchy problem. It 
also relates fermion masses to the CKM and MNS matrices, an appealing 
feature which has previously been considered in other contexts, e.g., [1]. 

By a rotating mass matrix we mean a rank-one mass matrix whose top 
eigenvector changes direction in generation space as a result of its evolution 
with renormalization scale. In this paper we wish to show that the rotating 
mass matrix hypothesis is fully consistent with the current experimental data 
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on both the fermion masses and mixing angles, except for the masses of the 
hghtest quarks for which the imphcation of the hypothesis is not clear. This 
work represents a serious advance from previous tests, which were either 
model-dependent [2], or limited to 2 generations by available data at the 
time |3]; and none of those could incorporate the CP phase in the CKM 
matrix as we now can, using a newly suggested mechanism [IJ. 

We find a smooth trajectory for the top eigenvector which results in a 
very good CKM matrix, with all the relevant parameters within experimental 
error, giving a strong CP- violating theta-angle of order unity. We then apply 
these results to Higgs decay and some flavour-violating decays for a range 
of Higgs masses. The Higgs decay presents interesting and testable features 
consistent with the predictions of [S]. The flavour- violating decays are found 
to be all within experimental bounds. 

2 Consistency of the Rotating Mass Matrix 
Hypothesis with Data 

2.1 Preliminaries 

The use of a rotating mass matrix to generate lower generation masses and 
nontrivial mixing was suggested in [2] and further described in [^. Here 
we need only to very briefly recall its essential features, and introduce some 
notations and formulae to be used in the phenomeno logical fit. 

Following Weinberg [7] we put the mass matrix in a hermitian form with- 
out 75. If, in addition, it is rank-one, then we can write 

m = rriTCta.^ (1) 

where rriT is the non-vanishing eigenvalue of m and a is its normalised eigen- 
vector. We shall assume that only ex., identical for the up- and down- type 
quarks, changes under renormalization and not the eigenvalue rriT, which 
depends only on the type of fermions. In this work we are concerned only 
with a phenomenological fit and so do not refer to any specific renomalization 
scheme. Models for this were treated in [21 16]. For ease of presentation we 
shall now give formulae for the up- type quarks, so that rriT = mu. The cases 
for the down-type quarks and leptons are similar. 

In renormalization theory it is usual to define particle masses at a scale 
equal to the mass itself. Similarly, in the current scheme we define the state 
vector of the top quark, v^, as the massive eigenvector at the scale of the top 
quark, v^ = ott, where we denote cx.{^ = rrii) by ckj. The eigenvalue of the 



mass matrix at this scale will then give the top quark mass, so rriu = rrit- We 
then run the scale to the mass of the charm quark. The particle state vectors 
Vi and Vc must be orthogonal, but the eigenvector etc need not be orthogonal 
to vj. We thus take the charm quark state vector to lie in the subspace 
orthogonal to v^ in the direction of etc, i-e., Vc || projy± (cKc). At /i = mc the 
charm quark state vector, v^, has in general a component in the direction of 
the massive eigenvector ckc- Thus the particle state acquires some mass via 
a 'leakage' from the massive state at /x = mc, even though it was massless at 
/i = mt- One can think of an effective ck*^^ G v^, ol^^ = projyx(Q;). This 0:*^^ 
will correspond to the eigenvector of a 2 x 2 rank-one matrix (since at this 
scale only the up and charm quarks can be produced), whose eigenvalue at 
/i = rric will be rric. We then follow a similar procedure to define v„. This 
procedure gives the masses of the up-type quarks to be 

mt = mu, 

rric = mu{\ac-Vc\f, 

mu = mu{\ctu■^ru\f. (2) 

Thus all the quark states get nonzero masses, as experimentally indicated, 
while the mass matrix itself remains rank-one with two zero eigenvalues at 
any scale. This important property of the rotating mass matrix, which may 
at first sight be surprising, comes about because of unitarity which requires 
that the mass matrix be truncated to remove those states below the threshold 
at which these states cease to be physical states. Hence in this 'leakage' 
mechanism physical masses are eigenvalues of the physical (truncated) mass 
matrix and not necessarily eigenvalues of the rotating mass matrix, and so can 
be nonzero while the rotating mass matrix retains its two zero eigenvalues 
at all scales. The fact that we should truncate the mass matrix at scale 
thresholds (as at /^ = rric above) is similar to what one has to do to the 
multi-channel 5'-matrix [H [9]. At scales below the highest channel the 5- 
matrix will have to be truncated from say n x ra to (n — 1) x (n — 1) in order 
to maintain unitarity, while the ^-matrix itself remains analytic at all scales. 
For further discussions on these points see, e.g., PH] . 

It should be noted that due to confinement in the strong sector the rela- 
tion between experimental measurements of the light quark masses and the 
definition we adopt above is unclear. For this reason we do not include the 
masses of the u and d quarks in this study. 

As the rotation of ck will be driven by some renormalization group equa- 
tion we shall be interested in finding a smooth trajectory. Although in specific 
models [21 [6] so far constructed all particle types share the same trajectory, 
this is not necessary in general. We find, however, that we can well accom- 



modate all particle types on the same trajectory. 

The state vectors of both the up- and down-type triads are determined 
at scales equal to the masses of the second generation, so that the above 
procedure determines all six quark state vectors. This in turn gives a matrix 
V 

v., = (vrivf-^) (3) 

where |v^'' own)^ ^^^ ^^^ state vectors so obtained. This matrix V would be 
the CKM matrix if there was no extra contribution from the theta-angle term 
in the QCD Lagrangian, the effects of which will now have to be considered. 

It has long been known [11] that if there is at least one massless quark then 
the CP violating term in the QCD Lagrangian can be absorbed by a chiral 
transformation on the massless quark field. The quark must be massless since 
otherwise this chiral transformation would make its mass parameter complex 
which would in turn lead to CP violating effects which are not seen [T2] . 

Now the rotating mass matrix scheme has, as noted above, the special 
property that the matrix can remain rank-one, i.e., with two zero eigenvalues, 
while giving nonzero masses to all the quarks. This means that one can rotate 
away any theta-angle term in the QCD action by a chiral transformation 
without the necessity of having a physical massless quark [10] . Further, as is 
shown in [1], the phase removed by the chiral transformation gets transmitted 
by the rotation to the CKM matrix and gives rise there to a Kobayashi- 
Makawa CP-violating phase, even if the matrix in ([3]) is real. If we set 
up a Darboux triad consisting of the radial vector cx.{n), the tangent vector 
t(/x) II cx{^) and the normal vector i/(/x) orthogonal to both, then it is argued 
in [1] that in order to preserve hermiticity of the mass matrix the chiral 
transformation should be effected in the normal direction v{fi). If we choose 
axes such that u is the third axis, then the chosen chiral transformation will 
give rise to a phase rotation of the left-handed fields 

/I \ 

Po = 1 , (4) 

V e-^^/2 / 

where 6 is the strong CP theta-angle. 

If we consider the up- type quarks for a moment, we have that r and u 
lie in the same plane as v^ and v^ at fi = rrit, and so the Darboux triad must 
be related to the state vector triad by a rotation in this plane. 



^u = I cosujj —sinuu . (5) 



1 











COSUu 


— sin Uu 





sin uu 


COSUu 



Here uu is the angle between r and Vg. Following |1] we find that the CKM 
matrix is given by 

K^KM = {n^^p^nu)v{n-^^p^^QD) (6) 

where V is the unitary matrix in ([3j). The mass formula ([2]) is unaltered by 
this chiral transformation. For given V , Uu and ur, the experimental value 
of the Jarlskog invariant can then be used to fix 6'. 

Assuming that the phases in ot do not change with scale is equivalent 
to assuming that all CP violating effects in the CKM matrix come from the 
strong theta-angle term in the manner set out above. Since the masses and 
mixing angles depend only on the inner products of ex with the state vectors 
the phases will cancel out in this case. This means we can take these vectors 
to be real without loss of generality. For simplicity we shall do so in our fit 
below and V then becomes an orthogonal matrix. 

We could consider the leptonic sector analogously. However, it is not 
known whether neutrino oscillations violate CP symmetry. It is also not 
known if neutrinos are Majorana particles, which would cause extra phases 
to enter the MNS matrix. As experiments so far have nothing to say about 
these points, and very little is known about neutrino masses, we do not 
consider the neutrinos in this work. 

2.2 Fitting ol{^i) to data 

The rotating mass matrix scheme has a as a fundamental object, and the 
state vectors are derived from this. To test this hypothesis we shall start from 
real orthonormal state vectors and use experimental data to find a consistent 
trajectory of a. We first choose the up type quarks to have state vectors: 

V. = (1, 0, 0)t, (7) 

vc = (0, 1, 0)t, (8) 

vt = (0, 0, l)t, (9) 

which we are free to do. Once we have chosen a V^ we can define the down 
type quark state vectors, v^ ,,;, = Vwu^c,t- There are no explicit physical 
constraints on the matrix V so we are free to choose any orthogonal matrix. 
For the top quark ol{^ = rrit) = v^, and similarly Q;(/i = m^) = v^. The 
leakage mechanism then fixes cx.c and cUs in terms of the quark masses: 



OLc = \frrQmt Vc + a/1 - mdrnt Vt, (10) 

CKs = ^/mslmi, v^ + a/1 - ms/mb v^. (11) 

The two vectors v^ and v^ define a plane. All that the mass ratio mu/rrit 
tells us about ck^ is the angle which it makes with this plane. We can thus 
restrict a^, and similarly a^, to lie somewhere on a line, parametrised by 

t e [0, 27r): 



cXu = \/mu/mt v„ + a/1 - niu/mt Vccos(t) + a/1 - niu/nit Vtsin(t), (12) 
"d = \/md/mb Vd + a/1 - rrid/mb v^ cos(t) + a/1 - md/rrib v^ sin(t), (13) 

where we will choose the signs of the square roots so that we obtain a right 
handed triad. The restrictions on cx.i for the charged leptons were found in 
an analogous way, replacing {u,c,t) with (e,/i, r). 

Before we give any results we will take a moment to discuss their pre- 
sentation. The state vectors and a are in R^ and of unit length so take 
values on the surface of a unit sphere. We will represent their positions by 
stereographically projecting onto R^. It turns out that a does not need to 
rotate very far from fi = rrit to fi = rUe so most of the action happens in a 
small area on the sphere. We have chosen the south pole of the projection to 
be at the position of v,-. This means that there will not be much distortion 
introduced by the stereographic projection; geodesies in this region on the 
sphere will be almost straight lines on the plane. Figure [1] shows the unit 
sphere with the region we will be interested in enclosed in a box. The curve 
within the box shows the best fit line we find and the point shows the south 
pole of the projection, v,-. The projection itself is shown on the right. The 
metric on the sphere is given by 

d^" = 7r—4—^2(du' + dv') (14) 

(1 + M^ + V^)^ 

for coordinates on the plane u and v. Over the boundary box in figure [T] 
the metric ranges from Aidv?' + dv'^) to 3.76((iM^ + df ^); there is a maximum 
distortion of a length in the stereographic projection of 3%. 

Though we assume that the trajectory of ol is universal, there is no phys- 
ical constraint on the relation between the quark and lepton sectors. We 
thus have the freedom to match these two sectors to give the smoothest tra- 
jectory. As mentioned previously we also have a freedom in choosing the 
matrix V . To find the best matching we ranged over, and then applied sim- 
plex optimisation at good regions in, the parameter space of quark-lepton 
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Figure 1: The state vectors and ex range over a unit sphere in [R^. These are 
conveniently represented as points on the plane under stereographic projec- 
tion. The point shown lies at the south pole of the projection. The box on 
the sphere is the boundary of the plot on the right, and the curve shows the 
best fit line found. 



sector matching matrices and quark mixing matrices, V . The positions of 
Ol■t,c,h,s,T,^l were found and projected onto the plane. We then fit a cubic line 
to these points on the plane using a non-linear least squares algorithm. Pre- 
vious work [3j has shown that the cumulative arc length between cx^ can be 
modelled by an exponential function at high scales. Accordingly, we then fit 
an exponential function to the cumulative arc length, excluding the strange 
quark and the electron. The strange quark was excluded since the interpre- 
tation of its intermediate mass is somewhat uncertain in this scheme. The 
electron was excluded firstly since we are only interested in the relationship 
at higher scales and secondly since we do not have good restrictions on the 
cumulative arc length, as oce can only be constrained to lie on a line. Fi- 
nally we use the Jarlskog invariant to fix the value of 9 and so determine 
the magnitudes of the elements of the CKM matrix. We chose the quark- 
lepton sector matching and the matrix V which produced absolute values of 
the CKM elements close to the experimentally measured magnitudes whilst 
maximising the R-Squared value of both the cubic and the exponential fit. 
The R-Squared value of a fit is given by 



R^ 



Ef hji - fif 
Ef ivi - yf 



(15) 



where Hi are the y coordinates of A^ data points, fi are the y coordinates of 



the best fit hne at the same x coordinates and y = j^ X^i Vi- 

Figure [2] shows the cubic best fit hne along with the positions of cx^. The 
experimental errors in the masses of the quarks lead to an uncertainty in the 
position of olc and Oig. The 1 a errors in the quark masses restrict otc and 
OLs to lie on the lines shown. The cubic best fit line is 

0.75 a: -4.83x2 -9.19x1 (16) 

The cumulative arc length between ckx's was found after mapping the 
cubic best fit line onto the sphere through an inverse stereographic projection. 

Figure [3] shows that the arc lengths between cxf, ck^, 0:^-, ck^ and ck^ are 
well fitted by the exponential curve 

0.104exp (-1.228 logio(/i)) -0.0061 (17) 

for /i in GeV. This is in good agreement with the results from the planar 
approximation found in [5]. The point otg does not fit on the exponential 
curve. This is expected, as the light mass of the s quark means that the 
definition of mass we use is not entirely correct for it. The arc length to CKg 
is estimated by the arc length to oie line. This is just an estimate since otf. 
may sit anywhere on this line. It is however clear that CKg is unlikely to sit on 
the exponential curve found. Various models which employ the rotating mass 
matrix mechanism suggest a tanh(/^) like behaviour, with fixed points for the 
rotation at /i = and /i = 00, so it is not surprising that the exponential fit 
matches the data well at higher scales but not at lower scales. For the later 
work on Higgs decay we only need to model the behaviour at high scales. 

For this trajectory we find Ujj = 0.09 radians and Uo = 0.25 radians. 
Fitting a Jarlskog invariant of J = 3.05 x 10^^ gives a strong CP angle of 
1.45 radians. These results are in line with estimates in P]. The absolute 
values of the CKM matrix obtained are: 
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which can be compared with the experimental values [T3] : 

0.97419 ± 0.00022 0.2257 ± 0.0010 0.00359 ± 0.00016 
0.2256 ± 0.0010 0.97334 ± 0.00023 0.0415+|];[^[J}? | . (19) 
0.00874+°:°°°^6 Q Q4Q7 ^ Q QQ^g o.999133+[;';S^ 
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Figure 2: The positions of ol at various scales determined partly by experi- 
mental constraints and partly by our choices as described in the text. The 
position of a^ is constrained to lie somewhere on the dotted line. The cubic 
best fit line is shown by the solid line. The top shows these on the sphere, 
left shows them on an ellipsoid with the axes stretched to match those of the 
stereographic projection while the right shows the stereographic projection 
with the positions of olx indicated. 
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Figure 3: The cumulative arc length along the best fit line measured from 
OLt is well approximated by an exponential curve for all but ocg and CKg. A 
great circle gives an arc length of 27r in these units. 
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Figure 4: An illustration of the angle ojd {<^u)- Here the state vector triad 
and the Darboux triad have their origins at 015 (ctt). The solid line shows 
t{^) II Q:(/i) and the dashed line shows v^ (vc). It should be remembered 
that the axes are not of equal scale and that this is a stereographic projection 
of the vectors so the angles cannot be directly read off. 
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We find the unitarity angles, defined and measured [13] as 

a = arg(-^) = (88+^)°, (20) 

VcdV:A 1 ^^., 



(3 = aTg^-^j = -sin-^(0.681± 0.025), (21) 

7 = arg(-^) = (77+i)°, (22) 

to be a = 88°, sin(2/3) = 0.691 and 7 = 70°. 

3 Higgs Decay Branching Ratios 

Now that we have found a trajectory of a which agrees well with experimental 
constraints we can use it to ask what the rotating mass matrix hypothesis 
may say about Higgs decay. Though we have so far ignored the lightest 
quarks note that the up and down triads are fixed at /x = rric and n = rris 
respectively. We can make predictions for the branching ratios of modes 
involving the up and down quarks since we shall see that they depend only 
on the quark state vectors and the trajectory of ct for /x > mH/2. We can 
similarly make predictions for modes involving the electron. Previously [S| 
the following Yukawa couplings have been suggested, e.g., for the up type 
quarks: 

Ayk = puo^i'L4>ipROL'^ + h.c. (23) 

Choosing the gauge in which 0, the Higgs doublet, is real and points in the up 
direction and expanding the remaining real component about its minimum 
value Cw, 'pR = Cw + H, we obtain the zeroth order mass matrix, equation 
([1]), with mu = PuCw, and the first order coupling matrix of the Higgs boson 
to the quarks as 

Y = puaal (24) 

Since this Yukawa coupling matrix depends on scale, /i, we need to take care 
in defining the details of Higgs decay. In constraining the trajectory of a 
we were considering fermion masses, whereas now we are considering Higgs 
decay. These are clearly related processes but in integrating the assumed 
underlying RGE the constant of integration may differ. It turns out that the 
correct calibration is a factor of 2 change in scale, i.e., ctH = cx.{fi = itlh/'^) 
[S] . We take the scale of Higgs decay to be that of the Higgs mass, as usual. 
We can now use the best fit line to find the Higgs state tensor, cxcx\ as a 
function of Higgs mass. 
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Figure 5: r{H — > xx)/T{H — )■ hh) for various final state particles as predicted 
by the standard model (left) and the rotating mass matrix hypothesis (right). 
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We then get the coupling for Higgs decaying into xy as 

A{H -^ xy) = pt\v^ ■ cxhW^v ■ olh\- (25) 

We can ignore any kinematic factors by considering ratios of branching ratios: 



(26) 



LEP found a lower limit on the Higgs mass of 114.4 GeV at 95% confidence 
level [13]. Various upper bounds have been given for the standard model 
Higgs boson mass. Here we plot up to niH = 260 GeV. 

Since now we have explicit formulae for the best fit trajectory of ck and the 
relation between arc length and scale, at high energies, it becomes a simple 
matter to find the ratios of branching ratios, given by fl26|) . for a range of 
Higgs masses. Figure [5] shows the standard model predictions for Higgs decay 
along with the predictions from the rotating mass matrix hypothesis. The 
standard model predictions were found using HDECAY [14J. We can see that 
the CC decay mode is heavily suppressed, in accordance with the estimate in 
[5] . This mode is suppressed since near p = 2mt the eigenvector a is almost 
orthogonal to Vc, reducing the branching ratio to cc. The suppression is to 
such a degree over the whole range of Higgs masses that detection should not 
require large statistics. The ss and /i"/!"*" modes are also suppressed, though 
to a smaller degree. 
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Figure 6: r{H — )■ xy)/T{H — )► hb) for various flavour violating decays as 
predicted by the rotating mass matrix hypothesis. Note that T{H — )■ yx) = 
r(if — )■ xy). Below around 220 GeV, indicated by the dotted lines, threshold 
effects will influence the tc and tu decay modes. 



Although we propose no dynamical mechanism the rotating mass matrix 
picture generically predicts flavour violating decays, shown in flgure El If the 
Higgs mass is large then tc and tu decay modes are possible. To be in line 
with the rest of the estimates, where the decays are well above threshold, we 
neglect threshold effects for these modes. Above 220 GeV we believe that 
these effects will be negligible and that the estimates shown for the branching 
ratios are realistic. The branching ratio for H — ?■ T~fi~^ is almost three orders 
of magnitude higher than that for H — )■ fi^fi^. This is in stark contrast to 
the standard model predictions where H — > r^yU"^ cannot occur at tree level 
and so has a small branching ratio. These modes may have a cleaner signal 
than H — )■ cc, making it easier to detect. Flavour violating effects at the 
levels predicted here have been shown in [5] to be consistent with existing 
experimental bounds on flavour violation. 

4 Conclusions 

We have found that the rotating mass matrix hypothesis can accommodate 
the CKM matrix with a CP violating phase and, with a few caveats, match 
the experimental data. The a^'s can be placed on a smooth trajectory so 
that they give the correct masses, for all but the light quarks, and give a 
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good CKM matrix. Fitting the Jarlskog invariant gives a theta- angle of 1.45 
radians. 

We have found that for the heavy quarks and the charged leptons a simple 
exponential fit can model the cumulative arc length for scales above ?7i^, 
though it is unlikely to be a realistic fit below this scale. This behaviour 
has been captured by both a phenomenological model (DSM [2]) and a field 
theory (MBSM [B]), which have rotating mass matrices and predict a tanh(/i) 
like behaviour with fixed points in the rotation at yU = and /i = oo. 

Since we do not yet know whether neutrino oscillations violate CP sym- 
metry, nor whether they are Majorana particles, the neutrinos place no con- 
straints on ct{fi) in the scheme as envisaged here. 

The cubic best fit line allows us to predict branching ratios for a range of 
Higgs masses. These give H ^ cc suppression, in line with a previous planar 
approximation [5], along with fi'fi'^ and ss suppression. We find a notable 
braching ratio to r"/!"*" and give limits for other flavour violating decays. 
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